Sturm-Schrodinger equations: formula for metric 

Miloslav Znojil 
Nuclear Physics Institute ASCR, 
250 68 Rez, 
Czech Republic* 

Hendrik B. Geyer 
Institute of Theoretical Physics, 
University of Stellenbosch and Stellenbosch Institute for Advanced Study, 

7600 Stellenbosch, 
South Africa'*' 

April 15, 2009 



*e-mail: znojil@ujf.cas.cz 
^e-mail: hbg@sun.ac.za 



1 



Abstract 

Sturm-Schrodinger equations Hip = E Wip with H ^ W and W ^ 7^ I 
are considered, with a weak point of the theory lying in the purely numerical 
matrix-inversion form of the double-series definition of the necessary metric 
operator 6 in the physical Hilbert space of states [M. Znojil, J. Phys. A: 
Math. Theor. 41 (2008) 215304]. This shortcoming is removed here via an 
amended, single-series definition of 0. 
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1 The origin of Sturm-Schrodinger bound- 
state equations 

It has already been known to Liouville pQ that after one pays due attention to 
boundary conditions an elementary change of variables may transform a given 
Sturm-Liouville eigenvalue problem 

Vi 1 ^) = E n ^\r) , #M = + V«(r) (1) 

into another Sturm-Liouville eigenvalue problem of just slightly more compli- 
cated form, viz., 

{x) = En W(x) (x) , =~ + V (2) (x) . (2) 

It is important to notice that both these equations are isospectral, i.e., we 
have E n = E^ = E n at all n = 0,1,.... Recently, we employed such 
an idea in ref. [2] (to be called paper I in what follows). We interpreted 
there eq. ([1]) as a standard Schrodinger equation of Quantum Mechanics, with 
iJW = T + where the differential operator T represented kinetic energy. 
In parallel, the change of variables (viz., r — > x and ip^ — > -i// 2 )) has been 
prescribed in advance so that we arrived at explicit formulae for the two new 
functions W(x) and V^ 2 '{x) entering the new differential-equation form ([2]) of 
our quantum bound-state problem. In such a context we found it very natural 
to call eq. (|2|) a Sturm-Schrodinger equation with Hamiltonian H and weight 
operator W. 

The principal advantage of our slightly counterintuitive transition from 
eq. ([T]) to eq. ([2]) [notice that we would have P^fr) = 1 in eq. ([T])] has 
been found in a perceivable simplification of the physical asymptotic boundary 
conditions. This implied that also the practical (say, numerical) solution of 
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eq. ([2]) proved much easier and quicker than the solution of its predecessor ([I]) 
(cf. refs. [2] or [3] for more details). 



2 The emergence of a nontrivial metric in 
Sturm-Schrodinger bound-state problems in 
Quantum Mechanics 

One must be aware that the presence of the nontrivial weight operator W(x) 
may make the solution of eq. ([2]) less routine. In paper I we paid particular 
attention to models and ^fy where in spite of the reality of all the spectrum 

( F) 

of energies a naive choice of the two alternative representations n)^ L of the 
Hilbert space of states of our quantum system in question happened to lead to 
inconsistencies like negative probabilities etc. For this reason the superscript 
( F ' stands here for "false". 

In paper I we emphasized that the superscript ^ could equally well stand 
here for "friendly". Indeed, once we temporarily forget about inner products 

( F) 

and consider just the underlying two vector spaces L we may treat them 
simply as spanned by the respective eigenfunctions or generalized eigenfunc- 

m 2) (1 2) 

tions *'{■) = <#^> of our respective Hamiltonians. 

In such a setting the main difficulty as addressed in paper I was represented 
by our specific assumption of the manifest non-Hermiticity of our Hamiltonian 
operators, 

^(1,2)^^(1,2^ in H (F) =n (F^ (3) 

Fortunately, being well aware of the inadequacy of the underlying choice of 
the specific and most elementary (often called Dirac's or Lebesgue's) inner 
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products 

{^\4>) = J dtA(r)il>*(r)<Kr), |^) , \</>) € (4) 

we immediately recalled the well known recipe as described, e.g., in refs. [1] or 
[5j. 

Let us briefly recollect the mathematical core of the recipe which lies in a 
redefinition of the operation of Hermitian conjugation. Technically this means 
that in the first step we decide to treat our Hilbert space 1iS F \ formally, as a 
vector space of kets \if) complemented by the "usual" definition of the dual 
space [V^] (marked by a prime) of bras {tj)\. In this setting the operation 
of Hermitian conjugation 

^general) 

is, in general, non-unique [5J. Still, people 
usually ignore this ambiguity and decide to work in full analogy with the linear 
algebra and finite-dimensional Hilbert spaces. Thus, they specify the dual 
vector space of linear functionals, i.e., bra vectors as quantities constructed 
from their ket-vector partners by the very specific Dirac's operation Ti Dtrac ) 
of transposition plus complex conjugation. 

In the second step we replace the above-mentioned Dirac's definition of 
conjugation in Hilbert space , VIZ., 

(iP\ := T^ Dirac ^ \if) in H {F) (5) 

by a non-Dirac definition of conjugation which amounts to an introduction 
of another Hilbert space 7i^ s \ This new, unitarily non-equivalent (i.e., inner 
product non-preserving) representation of Hilbert space is to be formed by the 
same vector space = and by a different "primed" vector space of 

linear functionals [V 1 ^]' / [V^l . 

On the level of notation we must graphically distinguish between the dual 
elements (if>\ G [V (F) ]' given by eq. © and, in our notation [2], elements 
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(("01 € [V^] denned by the different relation 

((rp\ = (0| e : =T inon - Dirac ^\ip) in n (s) , e^i 

or, in the coordinate representation, 

((0|-)= fdfi(r)r(r)e(r,-) . (6) 



{ EA 

This redefinition realizes the transition from the false Hilbert spaces 7^ 2) 
to their unitarily non-equivalent alternatives 7£/f 2) wnere the superscript ^ 
stands for "standard". In the coordinate represenation one speaks about the 
modified, double-integral definition of the inner product 

((0|0) = jd^x) Jdn(y) r(x)e(x,y)4>(y) (7) 

(S) 

applicable to any two wave functions and in H 12 j. 



3 Generalized Dyson-type mappings Q 

The key purpose of introducing the above-described nontrivial metrics = 

(S) 

6(1,2) nes i n making our Hamiltonians self-adjoint in 7~^ 12 )- For our further 
purposes we shall drop the subscripts and restrict our attention just to the 
Sturmian cases with subscripts 2- Next we introduce an invertible map £1 : 
V — > A connecting our vector space of functions (with ket elements [0) € V) 
with another, abstract vector space A composed of certain not yet specified 
spiked-ket elements, 

|V^=O|0), |0>}eV, \^€A. (8) 

Using such a notation we shall assume the solvability of the non- self adjoint 
doublet of the Sturm-Schrodinger equations composed of eq. ([2]) and of its dual 
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defined in the sense of space i.e., in the same vector space V. Thus, in 

an obvious shorthand notation we shall consider the following two equations, 

H | A) = A W | A) , H ] \ A')) = \'W ] \ A')) (9) 

By assumption, they define two families of elements of vector space V. After 
the application of mapping O one also obtains the corresponding families of 
elements of the abstract vector space A. Working here, exclusively and solely, 
with the standard and usual Dirac's conjugation we may introduce also the 
dual space A' and a related third, abstract Hilbert space HS P ^ where the 
superscript ( p > stands for "physical". 

These conventions lead to several immediate consequences. Firstly, without 
any danger of confusion, the elements of the dual space of linear functionals 
may simply be denoted by the spiked ket symbols, 

-< VI = e A'. (10) 

Secondly, the lower-case isospectral equivalent h = QHQ,^ 1 of our original 
non-Hermitian upper-case Hamiltonians H ^ H> as well as the parallel partner 
w = QWQ -1 of any original non-Hermitian specific "weight" operator W ^ 

may and will be assumed self-adjoint in the corresponding abstract physical 
Hilbert space n^ p \ In a way explained in ref. [S] the latter space is specified 
by its spiked- ket elements flSJ) and by its spiked-bra linear functionals (|10|) . 

All this means that we must require that 

fct = (fr 1 ) 1 " #t fit = h, = (fr 1 ) 1 & = w , 

or, after a trivial re-arrangement, 

ijt = QHQ' 1 , W ] = QWQ- 1 (11) 
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where we abbreviated O = fit fi. This shows how the concept of Hermiticity 
in physical Hilbert space remains equivalent to the apparent non- or 

pseudo-Hermiticity or cryptohermiticity in TL^ F \ In the language of ref. [5J, 
property (jlip of all our operators should be called quasi-Hermiticity in space 
Tv- S ^ where the metric is nontrivial. 

The pull-back of all the theory to the abstract Hilbert space TL^ clarifies 
that all our present considerations need not leave the domain covered by stan- 
dard textbooks [6]. In particular, due to the assumption of invertibility of the 
mappings fi we may replace our two original upper-case Sturmian problems 
Q by their common lower-case reincarnation 

h\ \y = \w\ \y (12) 

which is necessarily self-adjoint in Ti^ p \ Its simplicity facilitates the derivation 
of the Sturmian orthogonality relations 

~<\\w\\'>- = -<\\w\\>- -<5a,a' (13) 

and of the Sturmian completeness relations, 

J = y^|A^ - , -. - — ^A|u; (14) 
<\\w\\>- 

A 

as well as of the Sturmian spectral representation 

/i = y^w|A^ — —, — — 7— <X\w (15) 

^ 1 ^A | u; | A^ 

A 

of the lower-case Hamiltonian in THS P K 
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4 Biorthogonal bases and the double-series 
formula for the metric 



Our original Hilbert space Ti^ is, by assumption, so simple that it makes 
sense to transfer all the above formulae to this space. Thus, the insertion of 
definitions (JSJ) and (fTOj) in eq. (fT3|) yields the orthogonality relations in HS F > , 

( A | ft f wQ\X') = ( A | W | A' ) = ( A | W | A ) • 6 X ,\> . (16) 

Similarly, the appropriately adapted version of the completeness is obtained, 

f -El*> (A |eV|A) < A|9ly - < 17 > 

A 

Finally, the spectral decomposition of the Hamiltonian acquires the following 
form, 

H = E^W we X wm (Mew. (is) 

Moreover, the double-ket eigenstates | A )) of W may be understood as equal 
to the products 9| A) [21 E]. 

The key benefit of our return to the metric-independent and computation- 
friendly space 1iS F ' is that we may evaluate and set the matrix elements 
( A | G W | A ) in eq. (|16p equal to one via a suitable normalization of the basis 
[5]. This convention will perceivably simplify our formulae. The simplification 
would remain applicable whenever the Hermitian product W stays posi- 
tive definite (which is to be assumed). The first consequence of the resulting 
simplification of the formulae is that we may introduce the further, curly- ket- 
marked wave-function vectors and their curly-ketket partner vectors defined, 
respectively, as follows, 

\^} = w\ip), | v }} = | v }) • 



In terms of these new abbreviations we may stay working inside the most 
comfortable Hilbert space 1iS F ' and write down the following alternative and 
optional forms of the orthogonality conditions, 

(A|eW|A'} = {{A|A') = ((A|A / } = «5 AiA /. (19) 

Similarly we may derive the two alternative forms of the completeness relations, 

/ = Ei a >« a i = Ei a hai. (20) 

A A 

Finally, the menu of several alternative spectral-representation expansions can 
be replaced by the more or less unique, most compact expressions 

w = "£\\}{{\\, h = ^\^}H{M (2i) 

A A 

representing the weights and the Hamiltonian operators. 

At this moment one could decide to employ the method of ref. [2] and to 
insert spectral formulae (I21j) in the fundamental quasi-Hermiticity constraint 
or equation W = H yielding 

E|A}}A{A|0 = E©|A}A{{A|. (22) 

A A 
This relation would strongly suggest that the most natural Sturmian analogue 
of the well known single-series W = I formula should be sought via the double 
series ansatz 

= E |A}}M AiV {{A'|, M\ ; x' = ( A | 6 | A') . (23) 

A, A' 

This approach has been accepted in paper I and contributed to its discour- 
aging conclusions. Fortunately, a better strategy exists. Its core and main 
consequences will now be described in our last section 
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5 The single-series formula for = Q(H, W) 



We arrived at the heart of our present message. The identity \ip)) = \ip) 
inspires us to start from eq. (|20p and to obtain the following unexpected but 
sufficiently simple result 

A 

The impression of an apparent non-Hermiticity of this asymmetric formula is 
misleading and it is virtually trivial to verify that O = 0t as required. 

We should add that as long as the less economical, double-series formula 
(|23p for the metric is concerned, it might still find some marginal applications 
in some more complicated quantum toboggans [7] etc. Moreover, its use may 
always be combined with our present, single-series recipe. For example, in 
comparison with the recommendations formulated in ref. [2], the necessary 
nondiagonal matrix coefficients are now much more easily evaluated, 

M\ t \i = < A | e | A') = ((A | A'} = {{XIW-^X') = ((A \ W~ 1 \ A' } . (25) 

Moreover, the single-series expansion (I24p of the Sturmian metric could prove 
too compact when the analysis of some symmetries is concerned. 

Of course, the advantages of our present single-series formula will prevail in 
the majority of applications where, typically, people truncate the infinite series 
in order to obtain a reasonable approximation of observable quantities [8], 
etc. Moreover, even without immediate numerical applications, the existence 
and compact form of such a formula can definitely be seen as an important 
indication of the mathematical as well as physical consistency and tractability 
of the whole family of non-Hermitian, cryptohermitian Sturmian models. 
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6 Summary 



The standard areas of applicability of the general concept of Sturmians has 
been reviewed by Rotenberg [9j. The scope of this review ranges from the 
electron-hydrogen scattering to the interatomic charge-transfer collisions and 
to the atomic-physics related solutions of the Faddeev and Born-Oppenheimer 
dynamical and phenomenological equations. In the not too remote area of 
mathematical physics the transition from bound states [cf. eq. (pQ)] to Sturmi- 
ans [cf. eq. ([2])] found typical applications in perturbation theory [10]. 




Besides such a role aiming at an immediate determination of energy-dependent 
couplings the phenomenological use of Sturmians involves mathematically mo- 
tivated constructions of the so called quasi-exactly solvable quantum mod- 
els [11]. The latter models with W(r) ~ r M at positive (half)integer M 
re-attracted attention to many open questions like, e.g., the completeness of 
the sets of the Sturmian wave functions [12J. Last but not least, one should 
keep in mind that certain exact Sturmians emerging at W(r) = r N with any 
N = —1, 0, 1, 2, . . . proved helpful, in the context of classical physics, in con- 
nection with the resonant internal boundary layer problem [13] . 

Our preceding paper I revisited the problem of Sturmians and extended the 
domain of their applicability to the models characterized by a loss of manifest 
Hermiticity of the Hamiltonian. It has been established there that there exists 
a close connection between the concepts of Sturmians and of the so called 
crypto-Hermitian Hamiltonians with real spectra. We emphasized there that 
the survival of the reality of the energy spectrum can be expected to occur not 
only in the non-Sturmian regime with W = I but also in the presence of the 
nontrivial weight operators W ^ I. 

In our present continuation of paper I we felt inspired by the lasting interest 
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in the applications of Sturmians mediated by the problems ranging from the 
very standard computational physics [14J to the rather exotic tobogganic mod- 
els using a fairly abstract concept of complexified coordinates [21 [7]. In this 
setting we perceive paper I as an introduction to the subject which, naturally, 
suffered from an incomplete understanding of many subtleties. 

One of these subtleties has been clarified in our present paper. Our main 
conclusion is that the parallels between Schrodinger equations of the respective 
types ([T|) and ([2]) are much closer than expected in paper I. In particular, we 
feel pleased by the fact that in contrast to the very sceptical expectations 
expressed in paper I we succeeded here in finding the entirely general formula 
for the metric O and in showing that this operator of key relevance can be 
expressed as a single infinite sum over certain elementary projectors at W = I 
as well as at W ^ I, i.e., in both the non-Sturmian and Sturmian cases, 
respectively. 
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